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Abstract. Atmospheres of evolved AGB stars are heavily affected by pulsation, dust formation and mass loss, and they 
can become very extended. Time series of observed high-resolution spectra proved to be a useful tool to study atmospheric 
dynamics throughout the outer layers of these pulsating red giants. Originating at various depths, different molecular spectral 
lines observed in the near-infrared can be used to probe gas velocities there for different phases during the lightcycle. Dynamic 
model atmospheres are needed to represent the complicated structures of Mira variables properly. An important aspect which 
should be reproduced by the models is the variation of line profiles due to the influence of gas velocities. Based on a dynamic 
model, synthetic spectra (containing CO and CN lines) were calculated, using an LTE radiative transfer code that includes 
velocity effects. It is shown that profiles of lines that sample different depths qualitatively reproduce the behaviour expected 
from observations. 
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1. Introduction 

Stars of low to intermediate masses at a late stage of their 
evolution are characterised by high luminosities, large exten¬ 
sions, low effective temperatures and pronounced variability. 
After core helium burning has ceased, they reach the evolu¬ 
tionary stage of the Asymptotic Giant Branch (AGB). Towards 
the end of this stage, repeated and explosive helium burning in 
a shell sets in (thermal pulses, TP). Nucleo-synthesis products, 
especially carbon but also s-process elements, can be mixed up 
from the stellar interior by deep-reaching convection (dredge- 
up) and change the atmospheric chemical composition (e.g. 
Busso et al. ll999l l. Starting with an oxygen-rich (C/0<1) atmo¬ 
sphere at the onset of the AGB-evolution, the stars can become 
carbon-rich (C/0>1) during the TP-AGB phase. This results in 
a metamorphosis of the spectral appearance: the spectral type 
changes from M to C. 

Red giants become unstable to radial pulsation with of¬ 
ten large amplitudes and long periods of a few 10 to sev¬ 
eral 100 days. Among these so-called Long Period Variables 
(LPV), Mira stars show a well-pronounced and regular vari¬ 
ability. While their amplitudes reach several magnitudes in the 
visual, they decrease toward the infrared and bolometric light 
changes are of the order of 1™'*®. 
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Pulsations of the stellar interior have a strong effect on 
Mira atmospheres as they trigger shock fronts, which prop¬ 
agate outwards. This causes a levitation of the atmosphere, 
which consequently becomes very extended. The occurrence of 
shock fronts has been suspected from observations of emission 
lines around phases of maximum light (e.g. Richter et al. IMTTI 
I2003> and by the detection of line doubling found for molec¬ 
ular lines in the near infrared (e.g. Hinkle et al. I1982> . In the 
cool, dense environments of red giant atmospheres, molecules 
can form efficiently, resulting in the typical line-rich late-type 
spectra. In the post-shock regions, dust grains can condense 
from these molecules (e.g. Sedlmavr l 19941 MillarLike 
the molecular species, dust condensates also depend on the 
chemical composition of the atmosphere. For oxygen-rich ob¬ 
jects, silicates (amorphous as well as crystalline) and refrac¬ 
tory oxides (e.g. Mg-Al-oxides) are found, while for carbon- 
rich ones amorphous carbon and SiC can be observed (Molster 
& Waters EQQ3}. Radiation pressure on dust particles is be¬ 
lieved to be a major reason for the slow stellar wind from AGB 
stars leading to high mass loss rates (dust-driven winds; see e.g. 
Simis & Wotike l2004t . A recent general overview of the com¬ 
plex phenomena occuring in AGB atmospheres was given by 
Gustafsson & Hofner J2()()4l and references therein). 

High-resolution spectroscopy in the near infrared (NIR) 
proved to be a very powerful tool to study dynamic processes 
throughout the atmospheres of these stars (Hinkle et al. I1982t . 
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as discussed in detail in Sect.|2 The aim of this paper is to re¬ 
produce the diverse behaviours of several molecular lines by 
modelling synthetic line prohles from dynamic model atmo¬ 
spheres. 

2. Observed line profile variations in AGB stars 

2.1. Radial velocities from spectral features 

Spectral energy distributions of AGB stars peak in the NIR re¬ 
gion around 1-3/tm. The spectra are densely populated by nu¬ 
merous atomic and molecular absorption lines, requiring high 
resolution for detailed spectroscopic studies. Radial velocities 
(RV) derived from shifts in wavelength of spectral lines can 
provide clues to the velocites in the layers where these lines 
originate. 

Line shifts in Mira spectra have been explored observation- 
ally for a long time. The first investigations of these variations 
were made in the blue and the visual spectral range (e.g. Joy 
119541 Wallerstein 119751 revealing puzzling results, as outlined 
in Hinkle (fTWSt and Hinkle et al. J1982t . Radial velocities 
from emission lines, observed around maximum brightness, 
were consistent with the picture of a pulsating photosphere in 
which shockfronts propagate. However no large amplitude RV 
variations over the light cycle for absorption lines were found, 
as would be expected if the spectral variability is coupled to 
radial pulsations. Only a systematic red-shift to the center-of- 
mass radial velocity (CMRV) for almost all phases was derived 
(e.g. Fig. 4 in Hinkle & Barnes n~979bt . As these expected peri¬ 
odic RV variations were then found for other atomic and molec¬ 
ular absorption lines in the red spectral region (Sanford ll950> 
and later in the NIR (see below), it was suspected that sufh- 
cient optical depth in the visual prevents the blue-shifted com¬ 
ponents of absorption lines from being observable (Hinkle & 
Barnes [T979bl Hinkle & Barnbaum ll996> . 

Major advances in the studies of LPV atmospheric kine¬ 
matics resulted from the availability of Fourier transform spec¬ 
trometers. By applying high-resolution spectroscopy in the 
NIR, new insights could be gained. This spectral region is es¬ 
pecially well suited for Miras, not only because their fluxes 
peak at these wavelengths but also because variability is less 
pronounced there compared to the visual (hence spectroscopic 
monitoring over the light cycle is easier, even during minimum 
phases). 

The oxygen-rich Mira R Leo was hrst studied in detail us¬ 
ing spectra in the range of 1.6-2.5;um for different phases. The 
results were presented in a series of papers by Hinkle dWTa 
and Hinkle & Barnes lll979al[T97%b . They analyzed the be¬ 
haviour of different molecular (CO, OH and H 2 O) and atomic 
lines in this wavelength region and derived RVs. Hinkle et al. 
( 119821 in the following HHR82) then presented a very thorough 
study for^Cyg, a Mira variable of spectral type S, with a time 
series of several spectra at L6-2.5/im (providing good cover¬ 
age of the pulsational period) and a few at 4.6jum. From these 
observations it is known that different lines behave differently. 
Originating in various depths of the very extended AGB atmo¬ 
spheres, RVs derived from line shifts provide information about 
gas velocities in these layers for different phases. By combin¬ 


ing results from divers spectral regions, Hinkle et al. were able 
to give an overall interpretation of the complex dynamics going 
on in the atmospheres of these stars during a lightcycle. 

The variation of NIR lines was then explored further 
(larger samples of objects, different types of LPVs, different 
molecules, etc.) by several authors, e.g. Hinkle et al. ifTWli or 
Lebzelter et al. (11999112001t . Wallerstein ( I1985t showed that 
most of the layers sampled by NIR CO lines can be associated 
with RVs of selected lines in the optical region and attempted 
to derive a stratigraphy for x Cyg at different phases by com¬ 
bining optical and near infrared data. Recently, observations of 
spectral variations in the visual range have been obtained by 
Alvarez et al. ( 120001 1^0 1 al 1200 1 bt . From these they deduce a 
tomography of LPV atmospheres and present statistics on line¬ 
doubling among Miras. 

2.2. Molecular lines of CO in the NIR 

Although lines of several molecules and atoms were analysed, 
CO lines proved to be most useful for these kinds of inves¬ 
tigations for several reasons: (i) both carbon and oxygen are 
abundant elements and CO has a large dissociation energy, (ii) 
CO molecules are formed and can be observed in all types of 
chemistry (M, S, C stars), (iii) being formed already in deep 
photospheric layers, CO is abundant and stable over the whole 
atmosphere and large parts of the CSE (eventually it is photo- 
dissociated by the interstellar radiation field), (iv) it experi¬ 
ences no depletion into dust, (v) three sequences of vibration- 
rotation bands with sufficient line strength are observable in 
atmospheric windows. The achieved results (especially from 
HHR82) demonstrate that different vibration-rotation bands of 
CO are ideal probes of various regions within extended Mira 
atmosphere. As HHR82 state, CO bands are useful, because 
the observable lines span a wide range of excitation energy 
and strength, and line lists of reasonable quality exist (line po¬ 
sitions, oscillator strengths). Also, individual, unblended lines 
can be identified in the NIR, as well as points of the continuum. 

CO Av=3 absorption lines, found in the H-band at 
/i^L6jum, show blue/red-shifts coupled to the photometric 
variability and even line-doubling around visual maximum. 
The typical S-shaped, discontinuous RV curves (e.g. Fig. 12 
of HHR82) can be regarded as an indication of atmospheric 
motion driven by pulsation and accompanied by shockfronts 
running through the photosphere. This behaviour can also be 
found for high-excitation lines of CO Av=2 in the K-band at 
/l»2.3jum, but the low-excitation lines there show smaller vari¬ 
ations not necessarily related to the pulsation cycle. While the 
first originate in pulsating layers, the latter are formed further 
out in the dust-forming layers. CO Av=l lines observable in 
the M-band at A^4.6pm sample the outflow in the outer parts 
of the atmosphere and show roughly the same RV at any phase. 
Figurel^gives an overview of this scenario. 

In Fig. 1 of Lebzelter & Hinkle a composite of RV- 

curves derived from CO Av=3 lines of all Miras observed so far 
was put together, showing a rather uniform picture (shape, am¬ 
plitudes of about 20-30 km s“' around the CMRV) independent 
of spectral type, period or metallicity. Apart from this common 
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characteristic, line profiles can become highly asymmetric; 
their shapes and variations appear different in each observed 
star. This is even more pronounced around 4 /rm (Lebzelter et 

al.lMm. 

2.3. Carbon-rich stars 

This paper is dedicated to the exploration of C-rich atmo¬ 
spheres, as the formation and evolution of C-rich dust is much 
better understood and thus implemented in the models used 
here (cf. Sect. Eland Hofner et al. lM)^ . 

Most of the stars observed possess oxygen-rich atmo¬ 
spheres. Much less observational material exists for Miras 
of spectral type C. After the early publication by Sanford 
jl950t . carbon-rich AGB stars were later studied by Phillips & 
Freedman (I1969> . Hirai J1983> . Keady et al. J1988> . Barnbaum 
Ill992alll992bt and Barnbaum & Hinkle (I1995t . S Cep is the 
only C-rich Mira for which extensive time series IR spec¬ 
troscopy has been obtained so far. The analysis was presented 
by Hinkle & Barnbaum J 19961 in the following HB96). S Cep 
was therefore chosen as the reference object for our modelling. 
Spectral features of CN are prominent in visual and NIR spec¬ 
tra of C stars and can in addition to CO be used to infer velocity 
information. 

3. Modelling line profile variations 

If the behaviour of different molecular features can be repro¬ 
duced by synthetic spectra, this would be an indicator of the 
quality of the models used and a confirmation of the correct¬ 
ness of the ideas about the dynamic processes going on in the 
outer layers of Miras. On the other hand, the interpretation 
of complex multi-component line profiles, caused by the non¬ 
monotonic velocity fields throughout Mira atmospheres, will 
be aided by such computations. 

To model the line profiles discussed above, an atmospheric 
structure is needed for which detailed radiative transfer in¬ 
cluding velocities can be solved to obtain synthetic spectra. 
Modelling the cool and very extended Mira atmospheres re¬ 
mains a challenging problem as they have complex, tempo¬ 
rally variable structures. Hydrostatic model atmosphere codes 
were only able to reproduce observed spectra for AGB stars 
with mild pulsations and higher effective temperatures, e.g. 
Jprgensen et al. J2000t or Loidl et al. J2001t for C-rich stars and 
Aringer et al. <19971120021 for O-rich stars. However, this is not 
an adequate approach for red giants with more pronounced pul¬ 
sation and mass loss. It is even impossible to fit observations for 
these objects at various phases by different hydrostatic models. 
Dynamic models are in particular required to reproduce phe¬ 
nomena such as Doppler-shifted spectral lines. 

Hill & Willson (I1979t created the first hydrodynamic at¬ 
mospheric model to describe the observed velocity changes 
and emission features. Since then, more sophisticated mod¬ 
els for the atmospheres of LPVs and dust-driven winds have 
been developed by different groups. Reviews on the different 
approaches were given by Hofner <19991 and Woitke <19981 
120031 . Considerable progress has been achieved in this field 


during the last few years through comparison with IR spec¬ 
troscopic data (ISO). In general, dynamic models are now ap¬ 
proaching quantitative agreement with observations (Hofner et 
al. l2003l Gautschy-Loidl et al. IMmI Hofner et al. |2004j. In 
contrast, there have been only a few attempts to model line pro¬ 
files and their variations for pulsating AGB atmospheres in the 
past, which shall be summarised in the following. 

3.1. Pulsating model atmospheres 

Based on simple synthetic atmospheric structures of Bessell et 
al. <1989t . Bessell et al. <1988> presented NIR CO lines that 
show asymmetries and line doubling as in the observations. 
Using the more advanced dynamic models of Bessell & Scholz 
<1989ll . Scholz <1992t showed how velocities in Mira photo¬ 
spheres distort line profiles of different species and affect mea¬ 
surements of equivalent widths and curves of growth for chem¬ 
ical analyses from these spectra. Starting with the same models, 
Bessell et al. <1996> for the first time calculated synthetic line 
profiles (CO Av=2, Fe I) for different phases, which showed 
qualitative similarities with some observed features (line-shifts, 
doubling, velocity amplitudes). Scholz & Wood <2000t then 
presented similar calculations for different lines (CO Av=2,3, 
OH Av=2) and derived conversion factors relating RVs derived 
from Doppler-shifted lines to actual gas velocities above and 
below the shock wave running through the atmosphere. 

A unique feature of the dynamical models by Bessell et 
al. <1996> and Hofmann et al. <1998> is that they are based 
on self-excited pulsation models and should therefore be suit¬ 
able to represent the inner pulsating layers of a Mira atmo¬ 
sphere. Bessell et al. (1996) fitted the temporal variations of 
sub-photospheric layers of the pulsation models with Fourier 
series and used these fits as inner boundary conditions for dy¬ 
namical calculations of the outermost layers with higher spatial 
resolution. Hofmann et al. <1998> directly used the structures 
from their new pulsation models but followed the same proce¬ 
dure as Bessell et al. for the older models. In both cases, the 
density and velocity structures of the dynamical models were 
used directly for the calculation of observable properties. Based 
on the densities, however, new values for the gas temperatures 
were derived using radiative transfer in about 70 frequencies 
before computing the detailed spectra. This was necessary be¬ 
cause the temperatures resulting from the grey dynamical cal¬ 
culations were too unrealistic. A drawback of this procedure 
is that it neglects the feedback of non-grey effects on the den¬ 
sity and dynamics which may be crucial. Furthermore, different 
sources of opacity data were used in the different steps which 
adds to possible inconsistencies. In contrast to that, the dy¬ 
namical models used in this paper are based on a simultaneous 
solution of hydrodynamics and frequency-dependent radiative 
transfer leading to consistent dynamical density-temperature 
structures. It should also be pointed out that our models in¬ 
clude a time-dependent description of dust condensation and 
allow for the formation of dust-driven stellar winds whereas 
the models of Bessell et al. ( fT9^ and Hofmann et al. ([T9^ 
are pure atmospheric models (no dust and therefore no mass 
loss). Also the latter are focused on O-rich chemistry; no dy- 

























































4 


Nowotny et al.: Atmospheric Dynamics in Carbon-rich Miras I. 


namic models for carbon stars have been computed by Bessell, 
Hofmann and collaborators so far. 

3.2. Dust-driven wind modeis 

Improving the semi-empirical approach of Keady et al. (IT^ 
for spectral synthesis. Winters et al. (I2000I I presented syn¬ 
thetic fundamental and first overtone CO line profiles and com¬ 
pared them with observed ones for the obscured, C-rich Mira 
IRCh- 10216. Calculating synthetic spectra for different phases 
of a set of dynamic models with different parameters, they 
demonstrated the influence of mass loss rates (and hence dust 
optical depth) on the shape of the profiles for observations at 
different pulsational periods. By suppressing CO absorption for 
different layers. Winters et al. investigated their contribution to 
the final profiles. 

The models used by Winters et al. ( I2000t are based on a 
detailed description of dust formation and stellar wind dynam¬ 
ics and are well suited for describing stars with heavy mass 
loss and spectral features coming from the optically thick and 
dusty outflow. This allowed them to interpret variations ob¬ 
served in the first overtone of CO as the formation and dynam¬ 
ics of discrete dust shells predicted by the models. On the other 
hand, these models contain only a crude description of the stel¬ 
lar atmosphere, using grey radiative transfer and no molecular 
opacities. This leads to considerable differences in the density- 
temperature structures, in particular in deep photospheric lay¬ 
ers and further out in the upper atmospheric layers where the 
dust formation starts. This may be one of the reasons behind the 
problems of fitting both the global spectral energy distribution 
and the line profiles with one consistent model, as mentioned 
by Winters et al. J2000t . They rescaled the density of the wind 
to get line profiles comparable with observations but used the 
original model to compute the spectral energy distribution. Due 
to the uncertatinties in mass loss rates derived from observa¬ 
tions, no detailed fitting of the observed SED with the purpose 
of constraining the densities was performed. Compared to the 
models of Winters et al., the dynamical models used here com¬ 
bine a more realistic non-grey description of the atmosphere 
with a similar treatment of dust formation and wind dynamics. 
This allows us to study also stars with moderate mass loss rates 
that are not completely obscured by dust, and lines which orig¬ 
inate in very different geometrical depths in the atmosphere or 
wind. 

4. Combined atmosphere and wind modeis 

The discussion in Sect.|3 shows that existing dynamical mod¬ 
els that were used to study line profile variations in AGB stars 
fall into two groups; atmospheric models (e.g. Scholz & Wood 
120001 or wind models (e.g. Winters et al. 12000b . The former 
deal with lines originating from the various layers of the pul¬ 
sating atmosphere, while the lines studied with the latter mod¬ 
els probe layers from the dust formation zone to the outflow 
region. The purpose of the new models presented in this paper 
is a consistent description of all these phenomena, i.e. simulta¬ 
neous modelling of lines originating in various layers (from the 
deep photosphere out to the dust formation region and beyond 


to the stellar wind region), with one single dynamical model 
for a given star. * 

For the calculations presented here, dynamic model atmo¬ 
spheres as described in Hofner et al. mm were used. Starting 
with a hydrostatic initial model, the equations of hydrodynam¬ 
ics, frequency-dependent radiative transfer and time-dependent 
dust formation are solved simultaneously to get an adequate de¬ 
scription of the highly dynamic AGB atmospheres. This results 
in a more realistic description of both the dust-free pulsating 
atmosphere and the dust-driven stellar wind compared to pre¬ 
vious models. The pulsation is simulated by a variable inner 
boundary below the stellar photosphere (piston). Luminosity 
changes sinusoidally due to this varying boundary condition. 
See Hofner et al. ( 12003b for more details about the modelling 
method. 

In contrast to O-rich objects, formation of polyatomic 
molecules and dust happens within similar temperature regimes 
for C-rich stars (Loidl et al. 11999b . It is therefore necessary 
that opacities of molecules (not included in the models de¬ 
scribed in Sect. l3.2b and dust (not included in the models de¬ 
scribed in Sect. l3.1b are simultaneously treated in the compu¬ 
tations of the model atmosphere to get realistic atmospheric 
structures and NIR spectra. In our current dynamical models, 
the formation, growth and evaporation of dust grains in the C- 
rich case is treated by the method described in Gauger et al. 
( 11990b . In the case of O-rich chemistry, on the other hand, dust 
is in our models at present simply included in the form of a 
parameterised opacity, due to uncertainties in connection with 
the mirco-physics of grain formation.^ Therefore we started our 
calculations with models for carbon-rich Miras. 

These models have proven already to represent observed 
properties of AGB stars rather well, such as the global shape 
of spectral energy distributions (Hofner et al. 12003b and low- 
resolution spectra over a wide wavelength range (Gautschy- 
Loidl et al. 12004b . Describing the temporal variability of var¬ 
ious molecular line profiles visible in high-resolution spec¬ 
tra with the same models is another crucial test of whether 
the models reproduce structures and dynamics of Mira atmo¬ 
spheres well. An important step is to reproduce simultaneously 
the observed behaviour of various types of lines (coming from 
regions with different velocities within the atmosphere) for sev¬ 
eral phases during the light cycle, and to compare synthetic 
RVs with observed ones. 

For the line profile modelling discussed in this paper, a 
representative model atmosphere for a typical C-rich Mira 
(like e.g. S Cep) was calculated. The chosen parameters are 
shown in Tabled The same model was also used for com¬ 
parison with observed low-resolution IR spectra of S Cep in 
Gautschy-Loidl et al. (12004b and can be found in their Table 1 
(110t26cl4u4f20pi). Figured shows the corresponding spatial 
structures; the radius coordinate is plotted in units of the stel¬ 
lar radius of the initial model (R,, calculated from L, and T,). 

' A common feature of all these dynamic models is that they are 
computed in spherical geometry; effects of asymmetry therefore can¬ 
not be studied with these ID-models. 

^ Self-consistent models for O-rich mass-losing Miras were pre¬ 
sented by Jeong et al. 
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Fig. 1. Radial structures of the initial hydrostatic model and 
three selected phases 0boi of the dynamic model atmosphere 
used here for calculating synthetic spectra. Plotted are gas tem¬ 
perature (a), gas density (b), degree of condensation of carbon 
into dust (c) and gas velocity (d). The convention of assign¬ 
ing negative velocities to outflowing matter was adopted from 
observational publications for compatibility (cf. Sect. 15.2> . 


While the hydrostatic initial model is relatively compact, the 
atmosphere becomes much more extended as pulsation leads 
to a periodic levitation of the outer layers (compare Fig. 2a 
of Hofner et al. | 2003 . Shockfronts triggered by the pulsa¬ 
tion propagate outwards. If temperatures drop below the dust 
condensation temperature and densities are high enough be- 



Fig.2. Movement of mass shells with time at different depths 
of the dynamic model atmosphere (cf. Hofner et al. 120031 
Fig. 2). Pulsation dominates the absolutely regular motion in 
the innermost parts («;0.8-1.8R»). Where infalling matter col¬ 
lides with outflowing gas from the next pulsation cycle, shock- 
fronts emerge and propagate outwards. Starting at ss2Rt, dust 
is formed in dense regions and radiation pressure on dust par¬ 
ticles triggers the stellar wind. In these regions velocities are 
no longer periodic. From «;4R* outwards the dust-driven wind 
dominates and the model shows a stationary outflow. 


hind the shockfronts, dust can form (note that the dust shells 
in Fig. B slightly lag behind the shocks in gas density seen in 
Fig .QJi). Radiation pressure acting on the dust particles leads to 
an even more extended atmosphere. The dust shells move out¬ 
wards, dragging along the gas by dynamic friction and a stellar 
wind develops. Figure[2 shows that the inner regions are domi¬ 
nated by the pulsation; the structures duplicate from one pulsa- 
tional period to the next. The situation changes, as dust is being 
formed, starting at around 2 R* (Fig. B). As the time scales of 
pulsation and dust formation are independent, the behaviour 
of these layers can differ greatly depending on the model pa¬ 
rameters chosen. Dust shells are formed and evolve in differ¬ 
ent intervals (dust-induced /c-mechanism, cf. Hofner et al. 119951 
or Fleischer et al. rr^ not necessarily periodically (compare 
Fig. 1 of Hofner & Dorfi I1997> and the atmospheric structure 
is no longer repeating for the same phases of succesive peri¬ 
ods. From r; 4-5 R, outwards, the model shows a steady out¬ 
flow with only small changes of the terminal velocity on time 
scales longer than the pulsational period. A transmitting outer 
boundary, fixed at 30R,, is applied. Figure|2lshows the tempo¬ 
ral evolution of different layers. 
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Atmospheric structures from three different, separated cy¬ 
cles - several snapshots from the temporal evolution of the dy¬ 
namic model - were used to calculate synthetic spectra. To de¬ 
note model structures and synthetic spectra, bolometric phases 
0boi within the lightcycle in luminosity will be used throughout 
the paper. 

Table 1. Parameters and resulting properties of the dynamic 
model atmosphere used in this paper, chosen to resemble the C- 
rich Mira S Cep. The notation was adopted from Hofner et al. 
<2003t : Awp - velocity amplitude of the piston at the inner 
boundary, {u) - mean outflow velocity at the outer boundary, 
{fc) - mean degree of condensation of carbon into dust. 


L, 

[Lo] 

lO"* 

M, 

[Mo] 

1.0 

T. 

[K] 

2600 

R. 

[Rq] 

493 


[AU] 

2.29 

Ig g. 


-0.94 

C/0 


1.4 

P 

[d] 

490 

Aup 

[kms^*] 

4.0 

AwJboi 

[mag] 

0.86 

(M) 

[Moyr-‘] 

4.3-10-'’ 

(u) 

[kms^*] 

15 

ifc) 


0.28 


5. Spectral synthesis - synthetic line profiles 

In the following section, the basic properties of synthetic high- 
resolution spectra based on the dynamic model atmosphere de¬ 
scribed above are presented. 

5.1. Spectral lines, opacities and radiative transfer 

For the line modelling, several molecular and atomic lines 
were chosen according to the following criteria: (i) used suc¬ 
cessfully in published observations, (ii) line lists for the spec¬ 
tral synthesis are available, (iii) lines are observable from 
ground (atmospheric transmission), (iv) line depths, (v) low 
contamination by other lines or pseudo-continuous opacity. 
The high-resolution spectral atlases of Hinkle et al. ([IMI 
[20001 . Wallace & Hinkle and Ridgway et al. (119841 

were useful for identifications. If necessary, in some cases these 
were also used to correct the positions of the molecular lines in 
the original line lists. 

For the spectral synthesis, opacities calculated from line 
lists were used for CO (Goorvitch & Chackerian I1994> and 
CN (Jprgensen & Larsson I1990t . as described in Aringer et 
al. (11997112002> . The spectra were computed with a resolu¬ 
tion of /l/A/l=300000 and then rebinned to 70000, which is 
comparable to observed FTS spectra. For the calculations we 
adopted Doppler profiles assuming a microturbulence velocity 
of 2.5 km s“', which is consistent with previous modelling (e.g. 
Aringer fl999t and should have no major impact on the compu¬ 


tations as it is small compared to the macroscopic gas velocities 
of the model. 

For consistency and in order to get realistic spectra, 
molecules other than CO and CN should be accounted for in 
the radiative transfer as well. Also the optical depth at which 
the lines originate should be as realistic as possible. Thus, the 
line list data for CO and CN was supplemented with opac¬ 
ity sampling (OS) data for other C-bearing molecules (as de¬ 
scribed in detail in Gautschy-Loidl et al. I2004II . given with 
much lower spectral resolution than that used in the compu¬ 
tation of the line profiles. In order to account for the quasi- 
continuous opacity contribution of these species, one opac¬ 
ity value^ - averaged from sufficient (»300) OS data points 
around the central wavelength of the considered line - was 
added to each wavelength point of the high-resolution spec¬ 
tra calculated with line lists. It should be emphasised that the 
method works only for molecules with many densely-spaced 
absorption lines affecting the spectra more or less continuously 
(’’pseudo-continuum”). As the absorption changes only slightly 
and a mean value is larger than the fluctuations, averaging de¬ 
livers meaningful numbers.'* We found that the only relevant 
molecule was C 2 H 2 . 

Furthermore, frequency-dependent opacities of C-rich dust 
(pure amorphous carbon, no SiC) were included using the data 
(set AC) of Rouleau & Martin ( I1991t . For the moderately vary¬ 
ing dust opacity, one mean, constant value was added at each 
wavelength point of the high-resolution spectra computed from 
line lists. 

The radiative transfer code which was adapted and used for 
the spectral synthesis is described in detail in Windsteig (I1998t : 
we only give a short overview here. The code solves radiative 
transfer in spherical geometry by integrating the frequency- 
dependent transfer equation along characteristic rays (a de¬ 
scription of the applied method can be found in Yorke fT988t . 
In the current version, velocity effects are also accounted for. 
Because of the diverse movements of atmospheric layers at dif¬ 
ferent depths, including the influence of velocities on the in¬ 
teraction between matter and radiation is by definition essen¬ 
tial to model the complex behaviour of spectral lines in Mira 
spectra. For all results presented here LTE conditions were as¬ 
sumed.^ Level populations were computed from Boltzmann 
distributions at the corresponding gas temperature (resulting 
from the non-grey radiative transfer of the atmospheric mod¬ 
elling) and then used to solve the observer’s frame radiation 
transport equation by means of a Rybicki scheme. The spatial 

^ Constant over wavelength for all layers, but variable as a function 
of depth. 

For other molecules (like C 2 ), the described method results in a 
large overestimation of the opacity, since opacity distributions with 
pronounced features must not be averaged. Opacities therefore cannot 
be approximated by the discussed pseudo-continuum approach. 

^ This assumption may not be justified for layers crossed by a shock 
wave and in the extended regions some R, away from the star. This has 
been explored e.g. by Schirrmacher et al. <2003> . who state that non- 
LTE effects "are generally more pronounced outside of the dust for¬ 
mation zone, where the densities are smaller”. It is supposed that these 
effects do not contribute to the studies presented here (also assumed 
by Winters et al. 120001 . 
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resolution of the radiative transfer (radial points and impact pa¬ 
rameters) is determined by the mesh points of the modelling 
and therefore temporally varying due to the adaptive grid used. 
Additional depth points are inserted in the case of large radial 
(velocity) gradients. First results obtained with this code based 
on dynamic models were presented by Windsteig et al. Ill998al 
ll998b|[T^ and Nowotny et al. Il2003ll2004t . 

All synthetic spectra are normalised relative to a computa¬ 
tion where only continuum opacity is taken into account (no 
contribution of the molecular line itself). 

5.2. Naming convention for veiocities 

In all observational papers, velocities are dehned from the ob¬ 
server’s point of view, which by convention is positive for ma¬ 
terial moving away from the observer and negative for mat¬ 
ter moving towards the observer. Thus blue-shifted absorp¬ 
tion lines - indicating outflow from the star - give negative 
RVs, while infalling matter revealed by red-shifted lines re¬ 
sults in positive RVs. For easier comparison, this convention 
was adopted for the plots here (e.g. Fig.[3. Line profiles are 
then plotted with RV as abscissa, computed from the central 
wavelength and by using the formula for Doppler shift. 

5.3. Infiuence of veiocities on opticai depth 

Figure|3demonstrates for the model used how velocities within 
Mira atmospheres can influence the process of line forma¬ 
tion, by considering the CO 2-0 R19 low-excitation line (cf. 
Sect. l5.5> as an example. In the lower panel dotted lines show 
how the values of radial optical depth increase inwards by cal¬ 
culating T from radial integration over opacities and densities 
for the rest wavelength in the line center (A) and two wave¬ 
length points in the line wings (B, C - see insert in the upper 
panel of the figure). Solid lines show the corresponding optical 
depths if the opacities in all layers of the model atmosphere are 
Doppler-shifted according to the gas velocity there. 

For the rest wavelength of the line (corresponding to 
Ry=0kms ' - A), a large fraction of the opacity is shifted 
away due to the outflow in the outer regions. Therefore t= 1 is 
reached much deeper in the atmosphere if velocities are con¬ 
sidered. It is the other way round for the two wavelength points 
in the line wings marked B and C. The locations where optical 
depths reach unity move further away from the star, because 
additional opacity is moved to these wavelengths by Doppler 
shift. The optical depth for B is almost constant between 1.5- 
2R,, as most of the opacity is red-shifted by strong infall. 
This redistribution of opacity results in a broadened line, with 
a blue-shifted line minimum. The strong influence of veloc¬ 
ities on the region of line formation is quite evident. If t= 1 
is considered to measure the approximate location of the line¬ 
forming region, this is shifted from R* [A] to ^2.5 R* [B] 
(corresponding to a difference of ^5.73 AU). 

The low-excitation CO Av=2 lines originate from atmo¬ 
spheric layers with a wide range of velocities, hence a broad 
profile, and large changes in t result. Other lines, e.g. CO 
Av=3, originate in a relatively narrow depth range with more 




0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0 

radius R [R.=493R0] 

Fig. 3. Upper panel: Gas velocities at a certain phase and syn¬ 
thetic CO Av=2 line profiles calculated with/without taking ve¬ 
locities into account in the radiative transfer. 

Lower panel: Demonstration of how velocities in the atmo¬ 
sphere change optical depths and the approximate region at 
which the line originates (indicated by the t= 1-line). 

uniform velocities. Therefore, the shifted and split components 
show roughly the same line width as profiles computed without 
velocities (Fig.|4}, and the difference of the position of t-1 is 
smaller if velocities are taken into account. 

5.4. Probing the puisationai iayers 
5.4.1. CO Av=3 lines 

Around \.6pm the continuous opacity has a minimum (H - 
peak), allowing us to see deep into the atmosphere. This, in 
combination with the large excitation energies of second over¬ 
tone vibration-rotation lines of CO at these depths, led to the 
fact that these lines have been most used to sample regions 
where motions are dominated by pulsation of the stellar interior 
(HSH84, Lebzelter & Hinkle El)02t . Figure 1 of HHR82 shows 
a clear demonstration of the very characteristic behaviour for 
the S-type Mira^Cyg by means of averaged line profiles for 
different phases. Coupled to the lightcycle, periodic variations 
in wavelength shift can be seen, which repeat in the same way 
every puisationai period. A blue-shifted component appears 
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Fig. 4. Time series of synthetic CO 5-2 P30 line profiles during 
one lightcycle with different spectral resolutions. 


around 0v~O.9, becomes stronger, moves towards red-shifts, 
becomes weaker and disappears approximately one period later 
at (^v~0.1. For a certain time interval around phases of visual 
maximum, the lines appear doubled. This is due to the fact 
that infalling material (red-shifted) can be observed together 
with gas being accelerated outwards by the next emerging 
shock wave (blue-shifted component) at the same time (com¬ 
pare Fig. 1 of Alvarez et al. I2000t . This behaviour results in 
S-shaped, discontinuous RV-curves as in Fig. 12 of HHR82 and 
is usually explained by radial pulsations and emerging shock- 
fronts. 

Since observing second overtone CO lines is difficult for 
C stars due to contamination of this spectral region by other 
molecules (mainly CN, C 2 ), no corresponding studies have yet 
been published. Nevertheless, we synthesized such line pro¬ 
files, since observational studies indicate a very similar dy¬ 
namical behaviour for spectral types M, S and C (compare 
Lebzelter & Hinkle l7()()2> . 

The CO 5-2 P30 line at 6033.8967cm-* (1.6573/tm) was 
chosen for modelling. CN and C 2 , which are also prominent 
in this spectral region, were not taken into account. Including 
opacities from line lists of these lines could (due to uncer¬ 


tainties in positions) influence the profiles by blending.® This 
specific synthetic CO line may then not be directly compara¬ 
ble to the corresponding observed (and by CN/C 2 influenced) 
one. But since all CO Av=3 lines in this region have similar 
excitation potentials and gf values, the synthetic line should 
be comparable to an observed average line profile (where the 
various influences smooth out). The synthetic spectra here are 
computed in a slightly different way than the preliminary re¬ 
sults presented and analysed in Nowotny et al. (I2004t . Dust 
opacities are from Rouleau & Martin J1991t instead of Maron 
( II99OI 1. The former contains more complete data over a larger 
spectral range and was also used for the modelling of the atmo¬ 
spheric structures. In addition, the pseudo-continuous opacity 
of C 2 H 2 is also included here, leading to smaller line depth due 
to a depressed ’’continuum”. 

Figure|4]shows time series of the synthetic second overtone 
CO line profiles for two different spectral resolutions. To get 
an idea of the typical line width, a spectrum computed for the 
phase of light maximum and without taking velocities into ac¬ 
count in the radiative transfer is plotted at the top. Due to the 
limited number of spectral points the profil looks asymmetric 
and not centered on the rest wavelength (RV-Q) at lower res¬ 
olution. In principle, the above-described typical observed be¬ 
haviour of CO Av=3 lines in M/S stars can be recognised from 
the synthetic spectra. 

Some deviations from observations of Mira variables can 
be found. Most noticeable is the fact that the transition from 
blue- to red-shift during phases (f)ho\~0.2-0.5 is only visible 
in the lower-resolution spectra. The line profiles look more 
complicated at higher spectral resolutions. There, the move¬ 
ment of the original component stops at RyaiOkms ' where 
it disappears, while another red-shifted component develops 
from (p\,o\~0.2 on. On the other hand, line splitting is more pro¬ 
nounced at higher resolution {(pboi-0.75). From the synthetic 
spectra it can be noticed that the (pseudo-) continuous opac¬ 
ity is strongest^ for 0boi~O.4, leading to a depressed ’’contin¬ 
uum” and apparently weaker lines. A more detailed analysis 
shows that the synthetic CO Av=3 lines emerge from depths 
of R=0.8-1.3Rt with gas temperatures of »2200-3500K. This 
seems consistent with the observational results of HHR82 
(their Fig. 5) or HSH84, listing excitation temperatures of 
ai2000^500K. As it would be expected for lines coming from 
the very inner parts of the atmosphere with periodic movements 
(Fig.|2}, line shapes are reproduced for all three chosen periods. 
Therefore only profiles from the first one are shown. 

5.4.2. CN Av=-2 lines 

From the visual range to beyond 2pm, spectra of C stars are 
dominated by features of CN and C 2 . Coming from deep pho- 

® Only from synthetic spectra of this spectral region, the CO 5-2 
P30 line appears to be in a ’’window” of the CN-line-forest. This may 
allow a comparison with observed high-resolution spectra in the fu¬ 
ture. 

^ Around minimum phase, the rate of dust production is highest and 
also the C 2 H 2 features are strongest for a dusty model like the one used 
here. 
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CN t\u = -2 synthetic 
R = 300.000 


CN A!/ = -2 synthetic 
R = 70.000 



Fig. 5. Evolution of synthetic line profiles for the CN 1-3 Q24.5 
line during a pulsational period at higher and lower resolution. 


tospheric layers (Fig. 3 of Loidl et al. ll999t . CN lines are also 
suited for investigating kinematics in pulsational layers, as it 
has been shown by HB96. HB96 used lines around 2.14//m, 
which are electronic transitions of the red system of CN with 
Av=-2. These lines show the same typical behaviour as found 
for the CO Av=3 lines (Sect. l5.4'ni . also leading to a discontin¬ 
uous RV-curve. 

For the spectral synthesis the 1-3 Q24.5 line at 
4871.340cm“' (2.0528/tm) was chosen. The influence of C 2 H 2 
is rather weak there, but observations are difficult due to telluric 
contamination. No molecules other than C 2 H 2 are relevant for 
the overall continuous opacity in this region. As in the case of 
CO, the behaviour of this CN line should be similar to those 
typically used in observational studies. 

Figurel^shows the results of the spectral synthesis for sev¬ 
eral phases during the lightcycle. As they are the same from 
one period to the next, spectra are again shown only for one 
lightcycle. Observed prohle variations can qualitatively be re¬ 
produced by the calculations. The pattern with blue-/red-shifts 
and line doubling around light maximum (Sect. lS.TJl l is dis¬ 
cernible. However, the line shapes appear somewhat complex. 
The splitting is less pronounced than for the CO Av=3 lines; 
the two components are visible but merge. The strength of the 
red component does not decrease continuously. Several com¬ 
ponents are visible at higher resolution (e.g. (phoX-O.SO). This 
is better seen by rebinning the synthetic spectra down to reso¬ 


lutions comparable to observed FTS spectra (70000). The line 
profiles are smoothed and especially the transition from blue- 
to red-shift during phases of xO.2-0.8 appears more distinct. 
While for example at phase 0.80 the two visible red-shifted 
components merge into one broad feature, it becomes difficult 
to measure the weak blue-shifted component. The (pseudo-) 
continuous opacity (C 2 H 2 , dust) is again strongest for (p\,oi~0.4. 

With line profiles being even more complex, CN lines qual¬ 
itatively show the same behaviour as CO Av=3 lines. But - still 
sampling deep layers driven by pulsation - they originate from 
slightly smaller optical depths. This can be deduced from line 
doubling at later phases. 

In addition, a CN Av=-1 line of the 0-1 red system in the 
H-band was calculated. This 0-1 Qi76.5 line at 6107.012cm ' 
(1.6374//m) shows almost the same behaviour (prohles, RVs) 
as described above, being only different in a very small shift 
in 0boi- This was also found for a Ti line in the K-band where 
observed line doubling around light maximum (HB96) could 
clearly be reproduced. 

5 . 5 . Probing the dust-forming region - CO Av=2: 

The first overtone vibration-rotation lines of CO in the K-band 
around 2fim are also clearly observable. These have been stud¬ 
ied in time series of spectra for several stars: for the C-rich Mira 
S Cep by HB96, for the C-rich Mira IRCh- 10216 by Winters et 
al. ( I2000> . for the S-type Mira ^Cyg by HHR82, for the O- 
rich Mira RLeo by Hinkle JWTSl l and other M stars e.g. by 
Lebzelter et al. 

For the line modelling a few suitable lines used in these 
studies were chosen, especially the CO 2-0 R19 low-excitation 
line at 4322.0657cm-' (2.3137/zm) and the CO 2-0 R82 high- 
excitation line at 4321.2240cm"' (2.3142/zm). Concerning the 
treatment of other molecules, the statements of Sect. EO are 
also valid here. Figure^ shows a few synthetic spectra based 
on our dynamic model. 

The behaviour of CO Av=2 lines is two-fold. The weaker 
high-excitation lines (e.g. R82 in Fig. |6|l act in the same way 
as second overtone lines and almost duplicate their S-shaped 
RV-curve. While this is clearly visible for M- (Hinkle I1978> 
or S-type (HHR82) stars, contamination by many other lines 
(mainly CN) in this spectral region hampers the line identifica¬ 
tion for C-rich stars (see HB96). Figure|^shows that our model 
can reproduce the scenario for high-excitation lines reasonably 
well (compare Fig.l of Alvarez et al. 120001 . 

In the following, only the stronger low-excitation lines (e.g. 
R19 in Fig.|^ will be considered. Easier to identify in the spec¬ 
tra, they show a somewhat different behaviour. Variability is 
less pronounced than for second overtone CO lines (which 
appear shifted/split with approximately the same line width). 
They show complex, asymmetric shapes and seem to consist 
of several components, which are not further separable (Fig. 3 
in HB96). This is reproduced by our model; Fig.0presents a 
series of line prohles. A comparison with the prohle calculated 
without taking velocities into account suggests that various lay¬ 
ers with different velocities contribute to the hnal, broadened 
shapes. The variation is not repeating for the same phase of 
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R19 R82 R18 



A [/Lirn] 

Fig. 6. Synthetic spectra (resolution of R-3QQ 000) containing 
CO Av=2 lines and demonstrating the different behaviour of 
low (e.g. 2-0 R19) and high excitation lines (e.g. 2-0 R82). See 
text for details. Marked are the rest wavelengths of the lines. 


different lightcycles (see Sect. Hi. Again the pseudo-continuous 
opacity (C 2 H 2 , dust) is strongest for (p^ox^OA. 

Observed radial velocities being more or less constant 
around the CMRV leads to the scenario of a static layer 
within Mira atmospheres (e.g. HHR82 or Tsuii FToMb - cf. 
Nowotny et al. il2(i0^ Paper II) for a more detailed discus¬ 
sion. Also for the synthetic line profiles in Fig.0 most of the 
time one blue-shifted main component around ^5 km s ' is vis¬ 
ible, which would be expected from the estimated region of 
line formation (Fig.|3- B) and the velocities there (Fig.l^. In 
these depths of »2-3 R,, gas temperatures of ss800-1500K are 
found. This seems consistent with the excitation temperatures 
of 800+100 K derived for the same lines in;^fCyg spectra by 
HHR82. Low-excitation first overtone CO lines can be used to 
probe dynamics in layers where dust is being formed and the 
outflow starts (Fig. Hd). Similar conclusions were also drawn 
by Winters et al. J2000t . 



CO Af= 2 low-exc. synthetic 
R=300.000 




-30 -15 0 15 30 -30 -15 0 15 30 

RV [km/s] RV [km/s] 



CO Ap= 2 low-exc. synthetic 
R=70.000 




-30 -15 0 15 30 -30 -15 0 15 30 

RV [km/s] RV [km/s] 


Fig. 7. Comparison of synthetic line profiles of the CO first 
overtone low-excitation line 2-0 R19 for comparable phases 
from different lightcycles. Substructures that are visible at 
higher spectral resolution {upper panel) are smoothed out at 
lower resolution {lower panel). 


5.6. Probing the outflow - CO Av=1: 

Finally we study the fundamental mode vibration-rotation lines 
of CO located at ^Apm. Not much observational material ex¬ 
ists, as observations in this spectral region are hampered by 
telluric absorption. Nevertheless, these lines are useful as they 
sample the outermost layers of Mira atmospheres. Results on 
this have been presented for^ Cyg by HHR82, for IRC +10216 
(CW Leo) by Keady et al. ( I1988t and for a Ori by Bernat et al. 
0929 ). 

The CO 1-0 R1 line at 2150.856cm-' (4.6493jum) was 
chosen for modelling. Only the pseudo-continuous opacity of 
C 2 H 2 and dust is taken into account, since no other molecules 
contribute significantly in this region. Figure]^ shows a time 
series of synthetic spectra compared to one calculated with¬ 
out velocities taken into account. The line profile looks qualita¬ 
tively the same for all phases of the model and shows a typical 
P Cygni-type shape with a deep, blue-shifted absorption com¬ 
ponent (from the outflowing material in the line of sight) and 
a superimposed red-shifted emission component (from the ex¬ 
tended regions around the star). The strength of the emission 


is variable; in Fig.|^the extreme cases are shown. Similar syn¬ 
thetic profiles were presented by Winters et al. J2000t and com¬ 
pared to the observed ones of Keady et al. {Mill. From radial 
gradients of t (as Fig.|3l it is found that these lines originate 
as expected in the wind region at »15R,. Low gas tempera¬ 
ture of »350-500K in the line forming regions there appear 
comparable with excitation temperatures of 300+200 K given 
by HHR82. 

6. Summary and outlook 

Extended AGB star atmospheres are complicated due to dif¬ 
ferent interacting physical processes such as pulsation, forma¬ 
tion of molecules and dust and eventually the developement of 
mass loss. This leads to complicated, non-monotonic veloci- 
tiy fields. From several studies in the past it is known that dif¬ 
ferent molecular lines in high-resolution IR spectra are useful 
for probing kinematics throughout these Mira atmospheres. To 
simulate and understand the processes taking place there, dy¬ 
namic model atmospheres are constructed. They are especially 
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CO Ai/= 1 synthetic 
R=300.000 


CO Lv= 1 synthetic 
R = 70.000 



Fig. 8. Synthetic CO fundamental lines 1-0 Rl, showing char¬ 
acteristic P Cygni-type profiles indicative of a stellar wind. 


needed to reproduce variable phenomena, one of the most cru¬ 
cial aspects being the temporal variation of line profiles. 

Using the model atmospheres of Hofner et al. (12003t we 
were able to reproduce the typical global velocity structure (de¬ 
rived from Mira observations) from the pulsating atmospheric 
layers through the dust-forming region out to the outflow re¬ 
gion (Fig.|3i. It was demonstrated that lines originating in lay¬ 
ers of different depths can be modelled simultaneously by only 
one consistently computed atmospheric model and successive 
radiative transfer. It is absolutely necessary to include veloc¬ 
ity effects in the radiative transfer to model the complex line 
profiles and their variations. 

The behaviour of lines sampling different regions, as known 
from observations, could qualitatively be reproduced, as sum¬ 
marised in Fig.|9] 

- CO Av=3 lines, CO Av=2 high-excitation lines and CN 
Av=-2 lines sample the deep photosphere, which is dom¬ 
inated by pulsation. At low spectral resolution (70000), 
time series of synthetic line profiles resemble observa¬ 
tions (compare Fig. 1 of HHR82), although some differ¬ 
ences remain (e.g. line doubling is significantly too weak). 
At highest resolution (300000), the profiles appear more 
complex. It will be interesting to see if future instruments 
(with higher resolutions than the FTS spectrograph used by 
Hinkle and collaborators) will allow us to verify whether 
there are such substructures in the profiles or if this is a 
pecularity of the model atmospheres used. 

- CO Av=2 low-excitation lines sample the dust-forming re¬ 
gion, where the stellar wind is triggered. They appear there¬ 
fore slightly blue-shifted (reminiscent of a static layer) and 
show complex, broadened shapes with irregular temporal 
variation. 

- CO Av= 1 lines probe the layers of steady outflow and show 
typical P Cygni-type shapes at any time. 



Fig. 9. Upper panel: Density structures of the dynamic model 
atmosphere for several phases during one pulsational cycle 
{(fthoX-O-O-l.0), demonstrating the larger extension compared to 
the hydrostatic initial model and the density variations - shock- 
fronts due to pulsation in the inner regions and dust shells prop¬ 
agating outwards further out. 

Lower panel: Corresponding velocity structures; approximate 
regions probed by different types of molecular lines are 
marked, inferred from the velocity behaviour and plots of the 
optical depths. 


Thus, it can be ascertained that the dynamic model atmo¬ 
sphere used here shows fundamental agreement with dynamic 
processes (pulsation, dust formation, mass loss mechanism) oc- 
curing in Mira atmospheres. 

A detailed comparison of the synthetic spectra with ob¬ 
served FTS spectra of the carbon-rich Mira S Cep and radial 
velocities derived from line profiles is presented in Nowotny et 
al. (1200^ PaperII). 

Acknowledgements. WN wishes to thank Kjell Eriksson and the 
Uppsala Observatory staff for their support (scientifically and finan¬ 
cially) during his stays at Uppsala. Sincere thanks are given to T. 
Lebzelter and J. Hron for careful reading and fruitful discussions. This 
work was supported by the “Ponds zur Fdrderung der Wissenschaftli- 
chen Forschung” under project number P14365-PHY and the Swedish 
Research Council. 





























12 


Nowotny et al.: Atmospheric Dynamics in Carbon-rich Miras I. 


References 

Alvarez, R., Jorissen, A., Plez, B., et al. 2000, A&A, 362, 655 
Alvarez, R., Jorissen, A., Plez, B., et al. 2001a, A&A, 379, 288 
Alvarez, R., Jorissen, A., Plez, B., et al. 2001b, A&A, 379, 305 
Aringer, B., Jprgensen, U.G., & Langhoff, S.R. 1997, A&A, 323, 202 
Aringer, B., Hofner, S., Wiedemann, G., et al. 1999, A&A, 342, 799 
Aringer, B., Kerschbaum, R, & Jprgensen, U.G. 2002, A&A, 395, 915 
Barnbaum, C. 1992a, ApJ, 385, 694 
Barnbaum, C. 1992b, AJ, 104, 1585 
Barnbaum, C., & Hinkle, K.H. 1995, AJ, 110, 805 
Bernat, A.P., Hall, D.N.B., Hinkle, K.H., & Ridgway, S.T. 1979, ApJ, 
233, L135 

Bessell, M.S., Brett, J.M., & Scholz, M., et al. 1988, in Atmospheric 
Diagnostics of Stellar Evolution: Chemical Pecularity, Mass Loss, 
and Explosion, lAU Coll. 108, ed. K. Nomoto, Springer, p.l87 
Bessell, M.S., Brett, J.M., Scholz, M., & Wood, PR. 1989, A&A, 213, 
209 

Bessell, M.S., & Scholz, M. 1989, in Evolution of Peculiar Red 
Giants, lAU Coll. 106, ed. H.R. Johnson, B. Zuckerman, 
Cambridge University Press, p.67 
Bessell, M.S., Scholz, M., & Wood, PR. 1996, A&A, 307, 481 
Busso, M., Gallino, R., & Wasserburg, G.J. 1999, ARA&A 37, 239 
Eleischer, A.J., Gauger, A., & Sedlmayr, E. 1995, A&A, 297, 543 
Gauger, A., Gail, H.-R, & Sedlmayr, E. 1990, A&A, 235, 345 
Gautschy-Loidl, R., Hofner, S., Jprgensen, U.G., & Hron, J. 2004, 
A&A, 422, 289 

Goorvitch, D., & Chackerian, C.Jr. 1994, ApJS, 91, 483 
Gustafsson, B., & Hofner, S. 2004, in Asymptotic Giant Branch Stars, 
ed. H.J. Habing, H. Olofsson, Springer, ch.4, p. 149 
Hill, S.J., & Willson, L.A. 1979, ApJ, 229, 1029 
Hinkle, K.H. 1978, ApJ, 220, 210 
Hinkle, K.H., & Barnes, T.G. 1979a, ApJ, 227, 923 
Hinkle, K.H., & Barnes, T.G. 1979b, ApJ, 234, 548 
Hinkle, K.H., & Barnbaum, C. 1996, AJ, 111, 913 (HB96) 

Hinkle, K.H., Hall, D.N.B., & Ridgway, S.T. 1982, ApJ, 252, 697 
(HHR82) 

Hinkle, K.H., Scharlach, W.W.G., & Hall, D.N.B. 1984, ApJ SuppL, 
56, 1 (HSH84) 

Hinkle, K., Wallace, L., & Livingston, W. 1995, “Infrared Atlas 
of the Arcturus Spectrum, 0.9-5.3 microns” (San Erancisco: 
Astronomical Society of the Pacific) 

Hinkle, K., Wallace, L., Valenti, J., & Harmer, D. 2000, “Visible and 
Near Infrared Atlas of the Arcturus Spectrum, 3727-9300 A” (San 
Francisco: Astronomical Society of the Pacific) 

Hirai, M. 1983, in The MK Process and Stellar Classification, ed. R. 

Garrison, David Dunlap Observatory, p.353 
Hofmann, K.-H., Scholz, M., & Wood, PR. 1998, A&A, 339, 846 
Hofner, S. 1999, in Asymptotic Giant Branch Stars, ed. T. Le Bertre, 
A. Lebre, C. Waelkens, lAU Symp. 191, p.l59 
Hofner, S., & Dorfi, E.A. 1997, A&A, 319, 648 
Hofner, S., Eeuchtinger, M.U., & Dorfi, E.A. 1995, A&A, 297, 815 
Hofner, S., Gautschy-Loidl, R., Aringer, B., & Jprgensen, U.G. 2003, 
A&A, 399, 589 

Hofner, S., Gautschy-Loidl, R., & Aringer, B., et al. 2004, in 
Proc. of ESO Workshop ”High-resolution IR spectroscopy in 
Astronomy”, ed. H.U. Kaufl, R. Siebenmorgen, A. Moorwood, 
ESO Astrophysics Symposia, Springer, p.279 
Jeong, K.S., Winters, J.M., Le Bertre, T., & Sedlmayr, E. 2003, A&A, 
407, 191 

Jprgensen, U.G., Hron, J., & Loidl, R. 2000, A&A, 356, 253 
Jprgensen, U.G., & Larsson, M. 1990, A&A, 238, 424 
Joy, A.H. 1954, ApJ SuppL, 1, 39 

Ready, J.J., Hall, D.N.B., & Ridgway, S.T. 1988, ApJ, 326, 832 


Lebzelter, T., & Hinkle, K.H. 2002, in Radial and Nonradial Pulsations 
as Probes of Stellar Physics, lAU Coll. 185, ed. C. Aerts, 
T.R. Bedding, J. Christensen-Dalsgaard, ASP Conf Series, 259, 
556 

Lebzelter, T., Hinkle, K.H., & Hron, J. 1999, A&A, 341, 224 
Lebzelter, T., Hinkle, K.H., & Aringer, B. 2001, A&A, 377, 617 
Loidl, R., Hofner, S., Jprgensen, U.G., & Aringer, B. 1999, A&A, 342, 
531 

Loidl, R., Lancon, A., & Jprgensen, U.G. 2001, A&A, 371, 1065 
Maron, N. 1990, Ap&SS, 172, 21 

Millar, T.J. 2004, in Asymptotic Giant Branch Stars, ed. H.J. Habing, 
H. Olofsson, Springer, ch.5, p.247 

Molster, F.J., & Waters, L.B.F.M. 2003, in Astromineralogy, ed. Th. 
Henning, Springer, p.l21 

Nowotny, W., Aringer, B., Gautschy-Loidl, R., et al. 2003, in 
Modelling of Stellar Atmospheres, ed. N.E. Piskunov, W.W. 
Weiss, D.F. Gray, lAU Symp. 210, F3 
Nowotny, W., Aringer, B., Hofner, S., et al. 2004, in Proc. of ESO 
Workshop ”High-resolution IR spectroscopy in Astronomy”, ed. 
H.U. Kaufl, R. Siebenmorgen, A. Moorwood, ESO Astrophysics 
Symposia, Springer, p.291 

Nowotny, W., Lebzelter, T., Hron, J., & Hofner, S. 2005, A&A ac¬ 
cepted (Paper II) 

Phillips, J.G., & Ereedman, R.S. 1969, PASP, 81, 521 
Richter, He., & Wood, PR 2001, A&A, 369, 1027 
Richter, He., Wood, P.R, Woitke, R, Bolick, U., & Sedlmayr, E. 2003, 
A&A, 400, 319 

Ridgway, S.T., Carbon, D.E., Hall, D.N.B., & Jewell, J. 1984, ApJS, 
54, 177 

Rouleau, F, & Martin, P.G. 1991, ApJ, 377, 526 
Sanford, R.F. 1950, ApJ, 111, 270 

Schirrmacher, V, Woitke, P, & Sedlmayr, E. 2003, A&A, 404, 267 
Scholz, M., 1992, A&A, 253, 203 
Scholz, M., & Wood, P.R. 2000, A&A, 362, 1065 
Sedlmayr, E. 1994, in Molecules in the Stellar Environment, lAU Coll. 
146, ed. U.G. Jprgensen, Lecture Notes in Physics 428, Springer, 
p.l63 

Simis, Y., & Woitke, P. 2004, in Asymptotic Giant Branch Stars, ed. 

H.J. Habing, H. Olofsson, Springer, ch.6, p.291 
Tsuji, T. 1988, A&A, 197, 185 
Wallace, L., & Hinkle, K. 1996, ApJS, 107, 312 
Wallerstein, G. 1975, ApJS, 29, 375 
Wallerstein, G. 1985, PASP, 97, 994 

Windsteig, W. 1998, Master thesis, Vienna University of Technology, 
Austria 

Windsteig, W., Hofner, S., Aringer, B., & Dorfi, E.A. 1998a, in Proc. 
of ESO Workshop ’’Cyclical Variability in Stellar Winds”, ed. L. 
Kaper, A.W. Eullerton, ESO Astrophysics Symposia, Springer, p. 
308 

Windsteig, W., Loidl, R., Hofner, S., & Dorfi, E.A. 1998b, in 
Fundamental Stellar Properties: The interaction between obser¬ 
vation and theory, ed. T.R. Bedding, lAU Symp. 189, p. 172 
Windsteig, W., Aringer, B., Lebzelter, T., & Hofner, S. 1999, in 
Asymptotic Giant Branch Stars, ed. T. Le Bertre, A. Lebre, C. 
Waelkens, lAU Symp. 191, P- 614 

Winters, J.M., Ready, J.J., Gauger, A., & Sada, P.V 2000, A&A, 359, 
651 

Woitke, P. 1998, in Cyclical Variability in Stellar Winds, ed. L. Kaper, 
A.W. Fullerton, ESO Astrophysics Symposia, Springer, p.278 
Woitke, P. 2003, in Modelling of Stellar Atmospheres, ed. N.E. 

Piskunov, W.W. Weiss, D.F. Gray, lAU Symp. 210, p.387 
Yorke, H.W. 1988, in Radiation in moving gaseous media, ed. Y. 
Chmielewski, T. Lanz, 18'* Advanced Course of the Swiss Society 
of Astrophysics and Astronomy (Saas-Fee Course), p.l93 



